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ABSTRACT 

We present new observational results obtained for the Galactic non-thermal radio source G328.4+0.2 
to determine both if this source is a pulsar wind nebula or supernova remnant, and in cither case, 
the physical properties of this source. Using X-ray data obtained by XMM, we confirm that the 
X-ray emission from this source is heavily absorbed and has a spectrum best fit by a power law 
model of photon index r = 2 with no evidence for a thermal component, the X-ray emission from 
G328.4+0.2 comes from a region significantly smaller than the radio emission, and that the X-ray 
and radio emission are significantly offset from each other. We also present the results of a new 
high resolution (7") 1.4 GHz image of G328.4+0.2 obtained using the Australia Telescope Compact 
Array, and a deep search for radio pulsations using the Parkes Radio Telescope. By comparing this 
1.4 GHz image with a similar resolution image at 4.8 GHz, we find that the radio emission has a flat 
spectrum (a w 0; S v oc v a ), though some areas of the eastern edge of G328.4+0.2 have a steeper radio 
spectral index of a ~ —0.3. Additionally, we searched without success for a central radio pulsar, and 
obtain a luminosity limit of £1400 <i$ 30mJykpc 2 , assuming a distance of 17kpc. In light of these 
observational results, we test if G328.4+0.2 is a pulsar wind nebula (PWN) or a large PWN inside 
a supernova remnant (SNR) using a simple hydrodynamic model for the evolution of a PWN inside 
a SNR. As a result of this analysis, we conclude that G328.4+0.2 is a young (< 10000 years old) 
pulsar wind nebula formed by a low magnetic field (< 10 12 G) neutron star born spinning rapidly 
(< 10 ms) expanding into an undetected SNR formed by an energetic (> 10 51 ergs), low ejecta mass 
(M e j < 5Mq) supernova explosion which occurred in a low density (n ~ 0.03 cm' 3 ) environment. 
If correct, the low magnetic field and fast initial spin period of this neutron star poses problems for 
models of magnetar formation which require fast initial periods. 

Subject headings: stars: neutron, stars: pulsars: general, ISM: supernova remnants, radio continuum: 
ISM, X-rays: individual 



1. introduction 

Stars with initial masses between ~9 and 25 M Q are 
expected to end their lives in a giant supernova (SN) ex- 
plosion during which neutron stars are created. The fast 
moving ejecta from the SN create a supernova remnant 
(SNR), while the particle wind produced by the neu- 
tron star as it loses rotational energy inflat es a pulsar 
wind nebula (PWN: [Gaensler fc Slang l2006ft . Initially, 
the PWN is inside the SNR - and when the PWN is 
detected i nside the SNR the system is called a "compos- 
ite" SNR (jHelfand & Beckedfl98l . The evolution of the 
central neutron star and the outer SNR affect the PWN, 
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and as a result the PWN goes through several evolution- 
ary phases while it is inside the SNR. This evolution is 
determined by the physical properties of the neutron star 
(specifically the initial period Pq, the braking index p , 
and the strength of the dipole component to the surface 
magnetic field B ns ), the SN explosion (explosion energy 
E sn and ejecta mass M eJ -), and the surrounding medium 
(ambient number density n). As a result, by measur- 
ing the properties of the PWN inside a SNR at a given 
time one is able to constrain these physical parameters 
which allows one to study the mechanisms behind both 
core-collapse SNe and massive star evolution. 
With this in mind, we present results of new radio 

1 The braking index is defined as f2 oc Q p , where f2 is the angular 
velocity of the neutron star's surface. 
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observations of this source with the Australia Telescope 
Compact Array (ATCA), as well as a new X-ray (X-ray 
Multi-mirror Mission; XMM) o bservation of G328.4+0.2 
(MSH 15-57; iMills et al.|[l96lh. This sour ce is a distant 
(d > 17 A ± 0.9 kpc: IGaensler et al.ll2000h . radio bright 
(flux density S u = 14.3 ± 0.1 Jy at v =1.4 GHz), po- 
larized, extended (diameter D ~ 5.0') radio source with 
a relativel y flat spectral index (a ~ —0.12 ± 0.03 where 
S v cx v a ] Gaensle r et al.l [2000h . Based on these radio 
properties, and the discovery of non-thermal X-ray emis- 
sion from this source by ASCA dHuehes et al.l[2000h . this 
source was classified as a PWN - the largest and most 
radio-luminous PWN in the Galaxy. In this interpre- 
tation, the expectation is that G328.4+0.2 is ~ 7000 
years old and powered by an extremely energetic neu- 
tron star (|Gaensler et al.ll2000f). However, follow up ra- 



dio polarimetry work ( Johnston et al.ll2004l ) implied that 
G328.4+0.2 is an older composite SNR in which the 
PWN is just a small fraction of the total volume, and 
as a result is powered by a significantly less ener getic 
neutron star then argued by IGaensler et all (|2000f ). In 
this paper, we analyze new observations of this source in 
order to determine the age of G328.4+0.2, the energet- 
ics of the neutron star and the progenitor SN, and the 
density of its environment. 

In S}2] we present new X-ray and radio observations of 
G328.4+0.2. In £j3j we first discuss the expected evolu- 
tionary sequence for PWNe inside SNRs ( 33. 1[) , making 
general comments regarding the expected observational 
signature of each phase. In §3.21 use the observational 
results presented in fJ3J to draw some initial conclusions 
about the nature of G328.4+0.2. In 2J we present a sim- 
ple hydrodynamical model for the evolution of a PWN 
inside a SNR, which we apply to G328.4+0.2 assuming it 
is a composite SNR ( N4.1.ip or a PWN ( 34.1.20 . Finally, 
in §3] we summarize our results. 

2. OBSERVATIONS 

In this Section, we present the data gathered in a 
XMM observation ( ffTTj) , a 1.4 GHz ATCA observation 
of G328.4+0.2 ( §2.2|) . and a search for a radio pulsar in 
this source ( §2.30 . 

2.1. X-ray Observations 

On 2003 March 9-10, G328.4+0.2 was observed for 
^50 ks by XMM. During this observation, the PN cam- 
era was operated in Small Window Mode, and the Mosl 
and Mos2 camera were operated in Full Frame Mode. 
The "Thick" optical filter was used due to the presence of 
numerous bright stars in the field-of-view of G328.4+0.2. 
The data were reduced with the software package XMM- 
SAS v 6.0.0 with calibration files current through XMM- 
CCF-REL-174, using the standard procedure for re- 
ducing XMM data outlined in the XMM-Newton ABC 
Guide 2 and the Birmingham XMM Guide 3 . 

2.1.1. Image Analysis 

A vignetting corrected 0.2-12 keV image from the 
Mosl and Mos2 instruments 4 , is shown in Fig. [TJ in 



http: / /hcasarc. gsfc.nasa.gov/docs / xmm /abc/ 



2 Available at 

3 Available at 

4 We did not use the PN data due to the substantially larger pixel 
size of this instrument. 



ttp: / /www. sr. bham.ac.uk/xmm2/guidc. html 



which we observe three spatial components to the X-ray 
emission: a bright, compact feature located along the 
SW edge of the X-ray emission ("Clump 1"), a fainter, 
slightly extended feature located NE of the compact fea- 
ture described above ( "Clump 2" ) , and extended diffuse 
emission, roughly 1' in diameter, surrounding the two 
features described above ("Diffuse"). From the "Clump 
1" region we detected 120 ± 12 counts above the back- 
ground between 0.5-10 keV in the Mosl detector and 
136 ± 12 in the Mos2 detector, from the "Clump 2" re- 
gion we detected 66 ± 9 counts in both the Mosl and 
Mos2 detectors, and from the "Diffuse" regions we de- 
tected 360 ± 20 and 380 ± 20 counts from the Mosl and 
Mos2 detectors, respectively. 

We determined the spatial properties of these com- 
ponents using the S herpa modeling software package 
(|Freeman et al.|[200ll) . Due to the low number of counts 
per pixel, we used the sim plex fitting method and min- 
imized the cash statistic (Cash 1979). We attempted 
to model Clump 1 and Clump 2 as circular 2D Gaus- 
sians, elliptical 2D Gaussians, or circular 2D Lorentzians, 
the latter of which is a good model for XMMs point 
spread function fPSF: iGhizzardi fc MoTendll 120021) . We 
attempted to model the Diffuse region as a circular or 
elliptical 2D Gaussian, and assumed a constant back- 
ground. We fit the observed image to all model com- 
binations of Clump 1, Clump 2, and Diffuse (attempts 
to eliminate one of these components resulted in signif- 
icantly worse fits), and the best fit was obtained for 
a model in which Clump 1 is a 2D Lorentzian while 
Clump 2 and Diffuse are elliptical 2D Gaussians. The 
fit parameters for this model are given in Table [TJ and 
the model and residuals are shown in Fig. [TJ and from 
this conclude that the emission from Clump 1 is consis- 
tent with the PSF of XMM. Additionally, from this fit 
we estimate that Clump 1 and Clump 2 are separated 
by ~ 10", while the centers of Clump 1 and Diffuse are 
separated by ~ 15". 

2.1.2. Spectral Results 

In generating the spectra, only events with flag = 
were used. Additionally, the event files were screened for 
background flares by binning the 10-15 keV light curve 
of each instrument by 50 s and then recursively flagging 
all bins with a count rate > 3cr above the average. This 
procedure removed 3.0 ks, 2.3 ks, and 0.6 ks of data from 
the Mosl, Mos2, and PN detectors, respectively. Spec- 
tra were extracted for the regions shown in Fig. [TJ and 
the resulting spectra were binned into a minimum of 25 
counts per channel and modeled using Xspec v 12. 2.0. 
The background regions used are also shown in Fig. [TJ 

To determine the composite spectra of G328.4+0.2, we 
jointly fit the spectrum obtained by the Mosl, Mos2, 
and PN detectors - shown in Fig. [2j The background- 
subtracted observed 0.5-10 keV count rate of G328.4+0.2 
was 0.012 ± 0.001 counts s _1 in the the Mosl detector 
(555 ± 29 counts), 0.012 ± 0.001 counts s" 1 in the Mos2 
detector (580±30 counts), and 0.045±0.001 counts s" 1 in 
the pn detector 5 (1576±63 counts). We fit the spectra to 
seven different models separately - a power-law, a black- 
body, bremsstrahlung, a Raymond-Smith plasma, and a 

5 The PN detector count rate does not account for 29% dead time 
since this instrument was operated in Small Window Mode. 
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power-law plus one of these three thermal models - all 
attenuated for interstellar absorption. Only the single- 
component models produced reasonable fits (reduced 
X 2 ~ 1), and the fitted parameters are presented in Table 
[1 Both the blackbody (kT BB ~ 1.7 keV, T BB ~ 20 MK) 
and the bremsstrahlung models (kT ~ 9 keV) require 
unrealistically high temperatu res, especially i f the X -rays 
are from a PWN as arg ued bv lHughes et all ((20001 The 
derived parameters for the power la w model are s imilar 
to that observed in other PWN (e.g. lGotthelJ[2003l ). and 
agree well with the results obtained for G328.4+0.2 by 
Hugh e~s~et al.1 (|2000f l. There is no evidence for thermal 
X-ray emission, which one would expect from a SNR, in 
this source. 

Since the individual regions discussed in £|2. 1 . II did not 
have enough counts for spectral fitting, we measured 
their hardness ratio (HR), defined as: 



HR = 



H-S 
H + S 



(1) 



where H is the number of counts in the Hard (higher 
energy) band and S is the number of counts in the Soft 
(lower energy) band, of the regions discussed in §2.1. ll to 
determine if there were any spatial variations in the X- 
ray spectrum of G328.4+0.2. Using the X-ray spectrum 
as a guide, we calculate HR with H as the number of 
counts between 4 and 8 keV and S as the number of 
counts between 2 and 4 keV. Since the pixels on the PN 
detector are sufficiently large that it is not possible to 
separate the emission from these regions, we only use 
data from the MOSl and MOS2 detectors. The calculated, 
background subtracted HR of G328.4+0.2 is 0.25 ±0.04, 
of the Clump f region is 0.29±0.07, of the Clump 2 region 
is 0.31 ± 0.11, and of the diffuse region is 0.24 ± 0.05 
(lcr errors). As a result, we conclude that there is no 
significant change in the X-ray spectrum of G328.4+0.2 
between these features. 

2.1.3. Timing Results 

A clear signature for the presence of a neutron star 
would be the detection of X-ray pulsations in the emis- 
sion from G32 8.4+0.2. We searched for this using the 
test defined bv lBuccheri et al.1 (|1983l ). where n is the har- 
monic number of periodic signal, for n = 1,2,3,4. The 
maximum frequency searched was t-Wx = l/(2nAt), the 
minimum frequency searched was 1/20 Hz, and the fre- 
quency step was l/4t D bs (5 x 10~ 6 Hz, oversampling the 
Nyquist rate by a factor of 2), where At is the time res- 
olution of the dataset (5.7 ms) and i bs is the length of 
the observation. We only used events from the PN in- 
strument (5.7 ms since it was operated in Small Win- 
dow Mode) due to the poor time resolution (2.6s) of 
the MOS data. Additionally, we only used events from 
the Clump 1 region because only only emission from the 
central NS should be pulsed and as the brightest X-ray 
region of G328.4+0.2, this region is the most probable 
location of any neutron star. Unfortunately, due to the 
large pixel size of the PN instrument this region is con- 
taminated by emission from Clump 2 and the Diffuse re- 
gion. The event times from the resultant event list were 
barycentered to the Solar System reference frame, and we 
searched for a periodicity over multiple energy ranges in 
order to maximize the sensitivity of our search. The most 
significant period detected was in the dataset which only 



included photons between 5 and 10 keV (159 photons), 
in which for n — 4 a signal with period P — 336 ms had 
Z\ = 49.5 in 2074786 independent trials for this value 
of n and energy range, which has a 12% chance of being 
a false positive, a < 2a result. Statistically, the most 
significant sinusoidal (n = 1) pulse was in the 1-20 keV 
dataset (340 photons), and had a period P = 73.1 ms 
with a Z\ = 34.4 in 8365396 trials, a 34% chance of 
being a false positive. Assuming that this signal is not 
significant, we derive an upper limit on the pulse fraction 
of 45% for a sinusoida l pulse profile and 22.5% for a 5- 
function pulse profile (jLeahv et al.lfl98 3i). Using the PN 
count rate and size of the Clump 2 and Diffuse regions, 
~ 1/3 of the counts in the Clump 1 region is contamina- 
tion from these regions. Accounting for this, we are only 
able to put an upper limit on the pulse fraction of 67% 
for a (5-function pulse profile. This upper limit is consis- 
tent with the pulse fraction observed from other young 
neutron stars. 

2.2. Australia Telescope Compact Array Observations 

Based on previous radio observations of G328.4+0.2, 
there was a dispute in the lite rature as to whether this 
source is a PWN, as argued bv lGaensler et al.l (12000D. o r 
a composite SNR, as argued by iJohnston et alT (2004). 
The argument for this source being a PWN centered 
on the flat spectrum of the radio emission, as well 
as the high degree of pola rization observed from the 
center (|Gaensler et al.ll2000D . while the radial polariza- 
tion angles observ ed at the edge of this is more consis- 
tent with a SNR ([Johnston et all l2004f ). If there is a 
SNR component in G328.4+0.2, we expect that some 
of the radio emission from this source should have a 
steep (a < —0.3) spectrum. To search for such emis- 
sion, we observed G328.4+0.2 for 12 hours at 1.4 GHz 
with the Australia Telescope Compact Array (ATCA) 
on 2005 June 25. Flux density calibration was carried 
out using an observation of PKS B 1934-638, and phase 
calibration was carried out with regular observations of 
PMN J1603-4904. The observation was carried out using 
two 128 MHz bands, one centered at 1.344 GHz and the 
other at 1.432 GHz, and the data reduction was done 
using the miriad software package. The observation 
was conducted when the ATCA was in the 6B configu- 
ration, which has a longest baseline of ~6000 m (~ 6'.'9) 
and a shortest baseline of ^200 m (~ 3.'4). As a re- 
sult, this dataset alone is not sensitive to large-scale 
emission from G328.4+0.2. To improve the sensitivity 
to diffuse emission, we combined this datas et with the 
1.4 GHz data used by iGaensler et al.1 (|2000[ ) as well as 
continu um data gathered in the Southe rn Galactic Plane 
Survey (jMcClure-Griffiths et aT1l2005[ K Total intensity 
images from this combined dataset were formed using 
natural weighting, multi-frequency synthesis, and maxi- 
mum entropy deconvolution. The final image, shown in 
Fig. EH has a resolution of 7'/0x5'.'8, and an rms noise of 
~ 0.15 mJy beam -1 . The measured 1.4 GHz flux density 
of G328.4+0. 2 is 13.8±0.4 Jy - con sistent with the value 
measured bv lGaensler et al.l (120001). F or a 4.5 GHz flux 
of 12.5 ± 0.2 Jv (]Gaensler et alj|2000h . this implies that 
G328.4+0.2 has a radio spectral index a = -0.03 ±0.03. 

As seen in Fig. [3J the radio emission from G328.4+0.2 
is very complicated, and contains multiple morphological 
features. The major features are: 
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• Central Bar: The Central Bar is the bright- 
est morphological feature in G328.4+0.2, and 
has been previously detected at both 4.8 GHz 
(iGaensler et al.ll2000D and 19 GHz (|Johnston et all 
I2004T ). The Central Bar runs roughly E-W, and the 
we stern edge of th e bar i s bifurcated, first noticed 
by Uohnston etldl (|2004h . The length of the bar 
is ~ 1.'75 long, and inside the bar there are three 
peaks in the radio emission. 

• Filamentary Structure A: These are the curved 
filaments near the center of G328.4+0.2 which ap- 
pear to be connected to the Central Bar, the bright- 
est of which is the "Y" shaped structure NE of the 
eastern edge of the Central Bar. In general, these 
filaments are more prominent on the eastern side 
of G328.4+0.2 and appear confined to the central 
region of G328.4+0.2, not extending much beyond 
the inner half. 

• Filamentary Structure B: These are the faint, 
radial filaments predominantly found on the west- 
ern side of G328.4+0.2, as shown in Fig. H The 
inner parts of these filaments are located ~ 2.'5 
from the center of G328.4+0.2, and their length 
varies across G328.4+0.2 - in the southern half, 
two filaments appear to extend to the edge of the 
source while in the northern and western parts of 
G328.4+0.2 they are substantially shorter. While 
some of these filaments are kinked or curved, most 
are fairly straight and radial in orientation. These 
feat ures were not detected b y lGaensler et all (|2000f ) 
and iJohnston et alJ (|2004h due to the insufficient 
u-v coverage of those observations. 

• Filamentary Structure C: As shown in Fig. [3l 
these are features located near the outer edge of 
G328.4+0.2 that are parallel to the outer edge. 
These features are more prominent and more plen- 
tiful in the eastern half of G328.4+0.2. 

• Outer Protrusions: The Outer Protrusions are 
faint features - the two most prominent of which 
are in the NE quadrant of G328.4+0.2 - that ex- 
tend beyond the outer boundary of G328.4+0.2. 
Several of these structures have bow-shock mor- 
phologies. 

The physical interpretation of these features will be pre- 
sented in H3.2I It is worthwhile to note here that several 
of these morphological features (e.g. the Central Bar and 
the internal filamentary structures) have been observed 
in other P WNe such as MS H 15-5 £ (jDickel et al.ll2000h 
and 3C58 (jSlane et al.ll2004f ). while others (e.g. the Fil- 
amentary Structure C and Protrusions) are more char- 
acter istic of SNRs, such as the Vela SNR (jBock et alJ 
[19981) . 

This XMM observation also allows, for the first time, 
a comparison between the radio and X-ray morphology 
of G328.4+0.2. As shown in Fig. [3j there is a signifi- 
cant offset between the X-ray and the center of the radio 
emission, with Clump 1 located ~ 80" from the center 
of the radio emission. Additionally, the extent of the 
X-ray emission is significantly smaller than that of the 
radio emission. A physical interpretation of the X-ray 



morphology and its relation to the radio emission will be 
discussed in £ 13.21 

2.2.1. Spectral Index Map 

The previous 1.4 GHz dataset had a resolution (~ 20") 
significantly worse than that of the 4.5 GHz data, and 
therefore is not suitable for determining if there are 
small scale changes in a inside G328.4+0.2. With these 
new, high-resolution 1.4 GHz observations, it was pos- 
sible to make a spectral index map of G328.4+0.2 us- 
ing our new 1.4 GHz data an d the 4.5 GHz data pre- 
sented by lGaensler et al.1 (|2000f ) since these datasets have 
comparable u-v coverage. To detect any variation in 
a, we made a spectral tomog raphy map of G328.4+0.2 
(jKatz-Stone fc Rudnickl 119971 ). To do this, we first 
produced a 1.4 GHz image of G328.4+0.2 from data 
matched in u-v coverage with the 4.5 GHz data, and then 
smoothed both the new 1.4 GHz image and the 4.5 GHz 
image to a resolution of 8" to account for the poorer reso- 
lution of the 1.4 GHz data. Finally, we produced a series 
of difference images (Idiff,a) using the following formula: 

/1.4\ Q 

Idi«,a = h.4 ~ 14.5 ( J^T J (2) 

where /1.4 and I4.5 are the 1.4 and 4.5 GHz images 
produced above. In this method, the spectral index of 
a region is determined by the spectral index at which 
it disappears from the difference image. As shown in 
Fig. [5l most of the radio emission from G328.4+0.2 has 
a spectral index between a ~ —0.1 and a ~ +0.1, while 
the outer edges of G328.4+0.2 have a steeper spectrum 
(a ~ —0.4) than the center, particularly the western edge 
of G328.4+0.2. This steeper spectrum material is coin- 
cident with some of the Filamentary Structure C dis- 
cussed in i j2.2[ but there are no spectral features associ- 
ated with any of the other radio morphological features 
in G328.4+0.2 or with the X-ray emission. A physical 
interpretation of these results will be discussed in ij3.2l 

2.3. Search for the radio pulsar at Parkes 

As part of a project to s earch for pulsar count erparts to 
all Galactic PWNe (e.g. ICamilo et all l2002aT ) . on 2005 
October 13 we observed G328.4+0.2 using the ATNF 
Parkes telescope in NSW, Australia. As for similar 
such work, we employed the central beam of the Parkes 
multibcam receiver at a central frequency of 1374 MHz, 
with 96 frequency channels across a total bandwidth of 
288 MHz in each of two polarizations. The integration 
time was 24 ks, during which total-power samples were 
recorded every 0.25 ms for off-line analysis. 

We analyzed the data wit h standard pulsar sea rching 
techniques using PRESTO (jRansom et all l2002t K We 
searched the dispersion measure range 0-2600 cm -3 pc 
(twice the maxi mum Galactic DM pred icted for this line 
of sight by the ICordes fc Laziol 120021 electron density 
model), while maintaining close to optimal time reso- 
lution. In our search we were sensitive to pulsars whose 
spin period could have changed moderately during the 
observation due to very large intrinsic spin-down. The 
se arch followed very cl osely that described in more detail 
in ICamilo et alJ (|2006f ). We did not identify any promis- 
ing pulsar candidate in this search. 

Applying the standard modification to the radiometer 
equation, for an assumed pulsation duty cycle of 10%, 
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and accounting for a sky temperature at this location 
of 15 K, we were nominally sensitive to long-period pul- 
sars having a period-averaged flux density at 1.4 GHz 
of 5*1400 > 0.05 mJy. In fact, this limit applies only 
to such long-period pulsars as to not be of practical 
interest for us: for a distance of ~ 17kpc along this 
line of sight, the expected DM is w 1200cm~ 3 pc, and 
the scatter-broadening of the radio pulses due to multi- 
path propagation is ex pected to be ~ 50 ms at 1.4 GHz 
(|Cordes fc Lazi o 2002). This would likely render pulsa- 
tions undetectable from any short-period pulsar, such as 
we expect to power G328.4+0.2, regardless of average 
radio flux. We therefore repeated the search at a higher 
radio frequency v, since the scattering timescale is ap- 
proximately cx i/~ 4 . 

On 2007 January 4 we observed G328.4+0.2 at Parkes 
at a central frequency of 3078 MHz, with 288 channels 
spanning a bandwidth of 864 MHz in each of two polar- 
izations. The total-power samples were recorded every 
1 ms for 30 ks, and analyzed in a manner analogous to 
that described previously for the 1.4 GHz data. This time 
a few somewhat-promising candidates were identified in 
the analysis, and a second 3 GHz observation was made, 
on 2007 March 19, for 36 ks. Analysis of this second ob- 
servation did not confirm the original candidates, and we 
have therefore not detected any radio pulsar counterpart 
for the PWN in G328.4+0.2. The sensitivity of our 3 GHz 
observations was about 0.03 mJy for long-period pulsars 
(P > 20 ms) and decreasing gradually for shorter peri- 
ods. For the predicted DM, at this frequency the scat- 
tering timescale is expected to be ~ 2 ms, comparable to 
the dispersion smearing across each individual channel, 
~ 1 ms. Propagation effects should therefore not have 
prevented the detection of signals with P > 5 ms. 

Converting the 3 GHz flux density limit to a frequency 
of 1.4 GHz, using a t ypical pulsar spectral index of -1.6 
(jLorimer et alJ fl995h . results in S1400 < 0.1 mJy. For 
a distance of ~ 17kpc, this corresponds to a pseudo- 
luminosity limit of L 1400 = 5i400d 2 < 30mJykpc 2 . This 
is comparable to L1400 of th e very young pulsars B1509- 
58, J1119-6127, and Crab (|Camilo et al.ll2002b| ). but a 
factor of about 60 greater th an for the young pulsar 
in 3C58 (jCamilo et al.l l2002d) . which has the smallest 
known radio luminosity among young pulsars. Based 
on these results, it is therefore entirely possible that 
G328.4+0.2 harbors an as-yet undetected young pulsar 
beaming toward the Earth with an ordinary radio lumi- 
nosity. 

3. INTERPRETATION OF X-RAY AND RADIO OBSERVATIONS 
OF G328.4+0.2 

The X-ra y spectrum of G 328.4+0.2 is characteristic of 
PWNe (see iGotthelJ 120031 for a compilation of the X- 
ray properties of PWNe), and therefore conclude that 
the X-ray emission from G328.4+0.2 comes from a PWN 
and not from a SNR. The same is true for the polarized, 
flat-spectrum radio emission detected from the center 
of G328.4+0.2 which are also characteristic of PWNe. 
Therefore, in the following discussion we assume that 
the X-ray and flat spectrum radio emission are both pro- 
duced by a PWN. In §3.1[ we discuss the evolutionary 
sequence of PWNe in SNRs and the observational signa- 
tures of each stage. In W3.21 we use the results from the 
observations presented in £12.11 and £12.21 to draw general 



conclusion about the properties of G328.4+0.2. 

3.1. Evolution of PWN in SNRs 

Both PWNe and SNRs are dynamic objects, and when 
the PWN is inside the SNR its evolution is affected by 
the behavior of both the central neutron star and the sur- 
rounding SNR. While the PWN is inside the S NR, it typ- 
ically goes through three evolutio nary phases (| Chevalier! 
1998; Ivan der Swaluw et al.l[2004h : 

• The Free- Expansion Phase — In this phase, the 
PWN freely expands into the cold material inside 
the SNR, sweeping up and s hocking the surround- 
ing e j ecta into a thin shell (IChevalier fc Franssonl 
H99llvan der Swaluw et al.ll2001h . Since the PWN 
is confined only by the shock wave its expansion 
drives into the surrounding SNR, and the velocity 
of this shock wave is much larger than the neutron 
star velocity, it is free to move inside the SNR with 
the neutron star. 

• Collision with the Reverse Shock As the 

SN sweeps up and shocks the surrounding a am- 
bient material, a reverse shock (RS) is driven 
into the ejecta. Eventually, the PWN will en- 
counter the RS, and as a result, can not con- 
tinue to freely expand inside the SNR because it 
is no longer in an essentially pressureless environ- 
ment. Initially, the pressure behind the RS is 
higher than the pressure inside the PWN, and as 
result the PWN is compressed. As it contracts, 
the pressure inside the PWN increases adiabati- 
cally and eventually will be higher than its sur- 
round ing, and will as a result re-expand inside the 
SNR (iBlondin et all 120011: iBucciantini et alJl2003t 
i Revnolds fc Chevalier! 119841 : Ivan der Swaluw et al.l 
120011 ). Once the PWN encounters the RS, the ex- 
pansion velocity of the PWN decreases signficantly 
and falls below that of the neutron star, which is 
unaffected by this collision. As a result, the neu- 
tron star can detach itself from its PWN. 

• Relic PWN Phase — When the neutron star 
detaches from the relic nebula, it forms a new 
PWN from the relativistic e + /e~ plasma it contin- 
ues t o inject into the SNR (|van der Swaluw et al.l 
I2004D . The PWN around the neutron star and 
relic nebula evolve differently. The relic nebula 
continues to contract/expand inside the SNR un- 
til it achieves pressure equilibrium with its sur- 
roundings, a process that can take many tens of 
thousands of years. The new PWN initially ex- 
pands sub-sonically, but when the neutron star is 
~ 2/3 of the way to the SNR shell, its velocity 
will become supersonic relative to the surrounding 
material and the PWN will take on a b ow-shock 
morphology (jvan der Swaluw et al.f 2004). Eventu- 
ally, the neutron star will leave the SNR, and as it 
passes through the SNR shell it may re-energize the 
surrounding SNR material, as possibly observ ed in 
SNRs G5.4-1.2 and CTB80 (jShull et al.lli989l ). 

As the PWN evolves inside the SNR, its appearance 
changes radically. During the Free-Expansion phase, the 
morphology of the PWN is determined by the properties 
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of the particle wind expelled by the neutron star. In 
these cases, the neutron star is in the center of the PWN, 
and there is no significant offset between the radio and 
X-ray emission from the PWN. In general, during this 
stage the PWN is located near the center of the observed 
SNR shell. Examples of PWNe in this evolutio nary stage 
are t he PWNe in SNR 0540-693 in the LMC (|Revnoldsl 
1 19851). as well as those i n the Milky Wa y SNRs G11.2- 
0.3 (iRoberts et al.ll2003t .iTam et al.l )2"00l) and G21. 5-0.9 
(|Matheson fc San-Harbll2005D . 

Due to the offset in the PWN's position with respect to 
the SNR's center as a result of the neutron star's veloc- 
ity or inhomogeneities in the ISM, it is expected that one 
side of the PWN wi l l encounter the RS before the othe r 
side ([Blondin et all 120011: Ivan der Swaluw et all |2004[ ). 
As a result, the PWN will no lo nger be symmetrically 
orien ted around the neutron star (|van der Swaluw et al.l 
l2004h . leading to an offset between the radio and X-ray 
emission from the PWN. Because the cooling time for X- 
ray producing electrons is very short compared to that of 
radio emitting electrons, the X-ray emission of a PWN 
is expected to be brightest at the current location of the 
neutron star while the radio emission reflects the effect 
of the RS on the PWN. The compression/re-expansion 
cycle triggered by the PWN/RS collision will also affect 
the appearance of the PWN. According to a spherically 
symmetric MHD simulation of a PWN in this phase, 
compression of the PWN by the RS leads to an over- 
pressurized region forming in the center of the PWN. 
Material injected by the neutron star after this point is 
then confined to the small region, leading to the forma- 
tion of a radio/infrared "hot spot" in the center of the 
PWN (|Bucciantini et al.ll2003f ). 

The PWN/RS interaction also leads to the forma- 
tion of hydrodynamic, primarily Rayleigh - Taylor (R-T) , 
insta bilities at the PWN/SNR interface (Blondin et al. 
2001). During the PWN's initial free -expansion phase, 
the shell of material swept up by the PWN is subject to 
both thin shell and R-T instabilities (|Bucciantini et all 
120041: Uunl 119981 ). but the growth rate of these features 
is expected to be sufficiently small that the PWN is not 
disrupted, especially if even a small percentage of the 
total energy of the pul sar's wind is in magnetic fields 
([Bucciantini et al.l [2004). However, during the PWN/RS 
interaction, rapid mixing of pulsar wind and SNR ma- 
terial is expected when th e PWN re-expands into the 
SNR ([Blondin et all l200lh . In fact, numerical simula- 
tions suggest that the PWN is disrupted after its first 
re-ex pansion into the SN R as a result of these instabil- 
ities (|Blondin et alj [20011 ). It is important to note that 
these instabilities are only expected to affect the relic 
nebula and not the new PWN formed by the neutron 
star further from the SNR's center. Once the composite 
SNR enters the Relic PWN phase, the X-ray emission is 
expected to be dominated by the new PWN since it con- 
tains the high energy particles recently injected by the 
neutron star, while the radio emission is dominated by 
the relic nebula which contains most of the older parti- 
cles that arc expected to contribute at low frequencies. 
As a result, during this phase the radio-emitting elec- 
trons are expected to be dominated by electrons injected 
during the free-expansion phase, while the X-ray is from 
electrons injected after the passage of the RS. 



3.2. Observational Results 

In this Section, we use the results from the observa- 
tions presented in 32. H and £12.21 as well as the basic evo- 
lutionary sequence for PWN in SNRs described above 
in H3.ll to make some initial statements on the nature, 
evolutionary state, and properties of G328.4+0.2. The 
discussion given below is a very general interpretation of 
observed radio and X-ray features in G328.4+0.2, it pro- 
vides a framework in which to test the various scenarios 
for G328.4+0.2 discussed in 34X11 and gXJ 

As mentioned in 3TJ there is a debate in the 
literature as to whether G328 4+0. 2 is a PWN 
(iGaensler et al. 2000; Hughes et al. 2000) or a composite 
SNR ([Johnston et al.ll2004h . In neither of the X-ray or 
radio observations presented above is there clear evidence 
(e.g. thermal X-ray emission or a bright, steep spec- 
trum, radio shell) for a SNR component. If G328.4+0.2 
is a composite SNR, then the outer boundary of the ra- 
dio emission likely marks the outer radius of the SNR 
component, while the observed flat-spectrum radio emis- 
sion and power law X-ray emission are emitted by the 
PWN. Using the extent of the flat-spectrum radio emis- 
sion shown in Fig. O to estimate the size of the PWN 
component, we obtain that the radius of the PWN in this 
source must be > 2/3i?G328, the radius of G328.4+0.2. 
If G328.4+0.2 is a composite SNR, then Filamentary 
Structure 3, which as mentioned in 32.2.11 might have 
a steeper spectral index than the rest of the radio emis- 
sion in G328.4+0.2, would be emission from the SNR. 
This emission is then possibly analogous to the corru- 
gated structures seen in . the NE part of the Tycho's SNR 
([Velazquez et alJ[l998t ). Additionally, in this case the 
Outer Protrusions mentioned in H2.2I maybe are ejecta 
"bullets", similar to th ose observed in t he Vela SNR 
(lAschenbach et al.lll995T ) and SNR N63A ([Warren et all 
I2QQ3J). 

If G328.4+0.2 is a PWN, then the outer boundary of 
the radio emission is the outer radius of the PWN and 
the SNR in which it resides is undetected, similar to the 
case for the Crab Nebula. As a result, the steep spectrum 
radio emission seen at the edge of G328.4+0.2, as well 
as Filamentary Structure C and the Outer Protrusions 
observed in the radio, correspond to material swept-up 
by the PWN. This is because radio emission from the 
pulsar wind is observed to have a flat spectrum, un- 
like this material. As a result, Filamenta ry Structure C, 
as wel l as the radial component seen by Uohnston et al.1 
(|2004h in the polarization angle along the outer edge of 
G328.4+0.2, are the result of the hydrodynamical (HD) 
instabilities at the PWN/SNR interface. Since these 
instabilities only occur when the PWN is accelerating 
the shell of swept-up material surrounding it, he current 
pressure inside the PWN (P pwn ) must be higher than 
that of the SNR material just of the PWN [P s „ r (i? pwn )]. 
Additionally, in this case the Outer Protrusions would be 
the result of the PWN currently expanding into a clumpy 
medium (e.g. the PWN analog of the process described 
for young SNR bv lJun et al.lfl996l) . which requires that 
the expansion speed of the PWN, v pwn , is currently posi- 
tive. In this case, it is possible that the SNR surrounding 
G328.4+0.2 will be detected at a later date, as was the 
case for G21. 5-0.9 ([Matheson fc Safi-Harbll2005T) . 

Regardless of whether G328.4+0.2 is a composite SNR 
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or a PWN, the offset between the radio and X-ray 
emission from the PWN component implies that the 
PWN/RS collision has already occurred. The Central 
Bar is then the remains of the over-pressurized region 
created when the PWN was compressed by the RS, with 
the bar-like shape of this region the result of either an 
asymmetric RS or the anisotropic wind emitted by the 
neutron star. The Central Bar should then consist of 
pulsar wind material, which accounts for the a ~ spec- 
tral index of this region as shown in Fig. [5J Additionally, 
the flat spectral index observed from Filamentary Struc- 
tures A implies that this feature is emitted by pulsar 
wind material. As mentioned in §3.1[ when the PWN 
re-expands after the initial compression by the RS, nu- 
merical simulations suggest that the rapid mixing of the 
SNR ejecta and pulsar wind materi al can then occur. In 
their 2D simulations of this process, iBlondin et all (|200lD 
observe features similar to that of Filamentary Structure 
A, and as a result we conclude that the existence of these 
features requires that the PWN has re-expanded at least 
once after its initial compression by the RS. Since the ra- 
dius of flat-spectrum radio emission from G328 is larger 
than the outer radius of Filamentary Structure A, the 
instabilities formed during this expansion must not have 
complet ely disrupted the PW N. This differs from the re- 
sults of IBlondin et all (|2001h . in which the PWN is dis- 
rupted during the first re-expansion. This is most likely 
due to the damping effect of the PWN's magnetic field 
on Ra leigh- Taylor instabili ties, which is not accounted 
for by IBlondin et all 1)200 If ). 

Finally, we discuss the X-ray emission seen from 
G328.4+0.2 which, based on its X-ray spectrum, is emit- 
ted by pulsar wind material. The observed offset between 
Clump 1, Clump 2, and the Diffuse regions of the X-ray 
emission implies that the PWN is not freely expanding, 
consistent with the explanation that the PWN has col- 
lided with the RS. We identify Clump 1 as the current 
location of the neutron star since it is the brightest X-ray 
feature, Clum p 2 as the location of the t ermination shock 
in the PWN (jKennel fc Coronitilfl98l . and the Diffuse 
emission is produced by recently injected plasma stream- 
ing away from the neutron star. In this case, the extent 
of the Diffuse component depends on the synchrotron 
lifetime of the X-ray emitting particles in the PWN. 

The X-ray emission also provides an estimate of the 
physical properties of the central neutron star, namely 
a measure of the neutron star's rotational spin-down en- 
ergy, E. A comparison of observed X-ray luminosity L x 
and E shows a trend that neutron stars with a higher L x 
have a higher E, and that the relation ship between these 
two quantities is (jPossenti et al .1120021 ): 

logL X)(2 _ 10 ) = 1.341og£- 15.34 (3) 
where £x.(2-io) is the X-ray luminosity of the source 
between 2 and 10 ke V, albeit with a significant scatter 
(jPossenti et al.ll2002| ). For the absorbed power-law fit to 
the X-ray emission from G328.4+0.2 given in Table 
the unabsorbed 2-10 keV flux of G328.4+0.2 is - 1 x 
10 -i2 ergg cm -2 g -i_ For a di stance to G328.4+0.2 of 

d = 17dn kpc, we obtain that: 

£x,(2-io) ~3.5d? 7 x 10 34 ergs s~\ (4) 
which, using Eq. ([3]), gives us an estimate of E: 

E ~ 1.7d} 7 49 x 10 37 ergs s _1 . (5) 



T his number is somew hat less than estimate obtained 
WGaensler et alJ (|2000T ) (E = 8.3 x 10 38 ergs s" 1 ), who 
assumed that E = 4 x 10~ 4 Lr, where Lr is the radio 
luminosity of G328.4+0.2. Our estimate of E is similar 
to the estimate bv lHughes et all (|2000Tl (E ~ 10 37 - 2 x 
10 38 ergs s _1 ), who also used the L x — E relation. 

4. SIMPLE HYDRO DYNAMIC MODEL FOR THE EVOLUTION OF 
A PWN INSIDE A SNR 

In order to determine what neutron star, SN, and am- 
bient density properties are required to produce a system 
with these properties described in i j3.2[ we have devel- 
oped a simple hydrodynamic (HD) model for the evo- 
lution of a PWN inside of a SNR, which we then ap- 
ply to G328.4+0.2 in g4~U T his model is b a sed la rgely 
on the models develo p ed by IBlondin et alJ (|2001h and 
van der Swaluw et all (|2001l ). The main goal of this 
model is to determine the radius of the PWN, i? pwn , 
as it progresses through the evolutionary sequence de- 
scribed in £|3 - 1 1 In this model, we assume that the PWN 
can be treated as a perfect gas with adiabatic index 
7 = 4/3 and that is expanding into a SNR filled with 
a perfect gas with adiabatic index 7 = 5/3. We also as- 
sume that the material swept-up by the PWN initially 
lies in a thin shell with inn er radius R — 23/24 R pwn 
(|van der Swaluw et al.l l200Th . as shown in Fig. [51 The 
dynamics of this mass shell are determined by the dif- 
ference in pressure between the PWN and SNR, and by 
calculating the radius of this mass shell we determine 
R P wn{t). Once the PWN enters the Relic PWN phase 
of its evolution, this model only determines the proper- 
ties of the relic nebula. What follows is a brief qualitative 
description of the model used, while the full suite of equa- 
tions used to implement it quantitatively can be found 
in Appendix lAl 

As mentioned above, we model R pW n(t) by calculating 
the outer radius of the mass shell swept up by the PWN, 
ignoring the effect of any instabilities which could dis- 
rupt this shell. We solve for R pwn (t) by assuming that 
we know the values for the relevant quantities at a time 
t — At, and then calculate them for a time t, since the 
relevant equations can not be solved at all times analyt- 
ically, we wrote a program in IDL to implement this 
numerically using the following procedure: 

1. Calculate R pwn (t + At) by assuming that the mass 
shell around the PWN between t and t + At moves 
with a constant velocity v pwn (t). 

2. Calculate the internal energy of the PWN, E pwn (t+ 
At), using the first law of thermodynamics: 

AE pwn —Et — Ppwn AVp W n (6) 

which takes into account energy losses from the adi- 
abatic expansion/contraction of the PWN as well 
as any energy input from neutron star into the 
PWN between t and t + At if the neutron star is 
still inside the PWN. Since we assume the PWN is 
filled with a 7 = 4/3 perfect gas, E pwn oc P pW n^pwn 
from the ideal gas law, where P pwn is the internal 
pressure of the PWN and Vp Wn is the volu me o f 
the PWN, as defined in Equations (fA5]and (|A"6j) . 
Since V pwn oc i? 3 , and 7 = 4/3 requires that 



mvu > P wn , we derive that P pwn oc i?~ 4 . As 
a result, if there is no input from the neutron star, 
then P pwn oc -Rp^n- The energy input from the 
neut ron star (AE psl ) is calculated by integrating 
Eq. (fAT)) between t and t + At. 

3. Calculate P pwn (t + At) using Equations (IA5I) and 

~ a. 



4. Calculate the pressure inside the SNR (P sn r), the 
density inside the SNR (p C j), the velocity of the 
material inside the SNR (u e j), and the sound speed 
of the material inside the SNR (c s ), at the outer 
radius of the PWN, R — R pwn (t + At) using a 
model for the evolution and structure of a SNR, as 
described in Appendix IA"1 

5. If the PWN is expanding faster than the SNR ma- 
terial around it, increase the mass of the shell sur- 
rounding the PWN, M SWi pwn(£ + At), accordingly, 
as described in Appendix |A"1 

6. Calculate the force on the mass shell surrounding 
the P WN, F pwn (t + At), using Eqs. (|AT4|) and 
(|A15[) . During the initial free-expansion of the 
PWN inside the SNR, these equation s redu ce to 
Equation A4 o f Ivan der Swaluw et al.1 (|2001h and 
Equation 14 of IChevalieri (|2005l) . 



7. Calculate the new velocity of the mass shell around 
the PWN, v pwn (t + At), assuming that any mass 
swept up by the PWN between t and t + At is done 
so inelastically: 

(t)+F pm {t + At)At 
v pwn[ t + m)- Msw (t + At) [n 

This is believed to be a reasonable approximation 
because the newly swept-up material is shocked by 
the mass shell and, as a result, its pre-existing mo- 
mentum is transferred to the internal energy of the 
mass shell. 

This model is assumes that both the SNR and PWN 
are spherically symmetric, the PWN remains centered 
on the center of the SNR at all times, the PWN has no 
effect on the evolution of the SNR, and that the ma- 
terial swept-up by the PWN is incompressible and has 
a negligible internal pressure. Additi onally, this model 
ignor es the effects of magnetic field (jBucciantini et al.l 
120031) and RT instabilities at the PWN/SNR in terface 
(|Blondin et~aT1l200H Ivan der Swaluw et al.ll200l on the 
properties of the PWN. Despite these simplifications, our 
model does a reasona bly good job o f repro ducing the re- 
sults for Model A in iBlondin et al.l (|2001l), as shown in 
Fig. [7} In general, relative to results of IBlondin et al.1 
(2001) and other authors, our model tends to result in 
larger oscillations in R pwn / R S m and a larger initial com- 
pression. The first discrepancy results from neglecting 
the effect of instabilities at the PWN/SNR interface that 
damp these oscillations, and the second from not includ- 
ing the effect of reflected shocks th at enter the PWN a t 
the time of the PWN/RS collision (jBlondin et al.ll2001l ). 
Additionally, we find that scenarios with the same to- 
tal amount of energy deposited by the neutron star into 
the PWN, P psr but with different neutron star properties 



(e.g. different values of Po and B ns ) produce the same 
behavi or of R uv/n (t). This i s different than the conclu- 
sion of Blond irTet alj ((2001), and believe that this dis- 
crepancy is the result of using a more realistic expression 
for E, Eg. (IA7I). tha n a step function, the form used by 
IBlondin et alj (|2001h . 

4.1. Application of Model to G328.4+0.2 

In the following discussion, we use the model given in 
2] for a PWN's evolution inside a SNR to examine the 
different possibilities for the nature of G328.4+0.2 given 
in We first anal yze th e possibility that G328.4+0.2 
is a composite SNR ( t|4.1.1|) . and then he possibility that 
G328.4+0.2 is a PWN (ST2j). This model requires six 
inputs: the characteristic timescale of the neutron star's 
spin-down, ro, initial spin-down power of the neutron 
star Eq, the velocity of the neutron star, v ns , the kinetic 
energy of the SN ejecta P sn , the mass of the SN ejecta 
M e j , and the number density of the surrounding material 
n. In order to calculate these values, we use the following 
information: 

• The distance to G328.4+0.2 is 17 kpc {d 17 = 1), 
which is the lower limi t on the distance to this 
source as determined by iGaensler et al.1 (|2000D us- 
ing Hi absorption. This implies that the current 
radius of G328.4+0.2 is P G32 8 = 12.5 pc. 

• The neutron star inside G328.4+0.2 is spinning 
down with a braking index p — 3, the braking in- 
dex produced by a pure dipole surface magnetic 
field. Additionally, we assume that the neutron 
star's moment of inertia is / = 10 45 g cm 2 , the 
value derived f or standard equations of st ate for 
a neutron star (|Shapiro fc Teukolskvl [l983h . Both 
of assumptions are standard in the literature (e.g. 
IBlondin et al.lf200lh . 

• The offset between the Clump 1, which as described 
in i j3.2l is believed to be the location of the neu- 
tron star powering the PWN, and the center of the 
radio emission is due to neutron star's spatial ve- 
locity, v ns . Using M2.1.H we determine that this 
observed offset corresponds to physical distance of 
~ 6.6di7 pc. Since the observed offset is due only 
to the neutron star's velocity in the plane of the 
sky, it is a lower limit on the true distance the neu- 
tron star has traveled since the SN explosion, r ns . 
If we assume that v ns — ^ns/^now? where t now is 
the age of G328.4+0.2, the observed offset allows 
us to estimate the minimum spatial velocity of the 



neutron star, 



equal to: 
6.6dn pc 

^now 



(8) 



Using Equation (|A19[) , we determine that for stan- 
dard initial periods (Po ~ 5 — 20 ms) and magnetic 
field strengths (P ns = 5 x 10 11 - 10 13 G), r varies 
from ~ 100 — 2000 years. To cover this range, we 
assume that To of the neutron star in G328.4+0.2 
can have one of three different values: 

r = 430, 770, and 1730 years. (9) 

which respectively correspond to a neutron star 
with B = 10 12 G and P = 5 ms, B = 3 x 10 12 G 
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and P a = 20 ms, or B = 5 x 10 11 G and P = 
5 ms. This range of m is similar to those used 
bv iBlondin et alj (l2001h, I van der Swaluw et al.l 



(20011), and lBucciantini et all (|200l) 



• G328.4+0.2 is expanding i nto a uniform m edium 
(s = in the notation of IChevalierl |1982j) . This 
assumes G328.4+0.2 is much larger than either 
the main sequence and late-stage wind bubble 
formed by its progenitor, both of which are ex- 
pected to have an interior p a r -2 density struc- 
ture. While the typical size of these structures 
is smaller than i?G328j this is not the the case 
for very massive stars (M > 15 Mq) for which 
these bubbles can reach sizes of ~ 100 pc or 
larger in low density (n < 1 cm" 3 ) environments 
(IChevalier fc Emmerinel 11989! Chevali er &: Liang] 
1989). In this case, our assumption of a constant 
density medium would not be correct. However, 
the effect of G328.4+0.2 still being inside a stellar 
wind bubble since this does not significantly mod- 
ify the evolution of the SNR. Since their no a priori 
information on the density around G328.4+0.2, we 
assume it is one of the assume it has one of the 
following values: 



logn 



-1.5,-1.0,0,0.5 cm- 



flO) 



which cover the range of densities in the warm ion- 
ized medium and the warm neutral medium. 

• The ejecta mass of the SN explosion that formed 
G328.4+0.2, M C j has one of the following values: 

M ej = l,5,lOM (11) 

and that the kinetic energy of the ejecta, E sn is: 

log(£ sn /10 51 ergs) = -0.5, 0, 0.5. (12) 

This range of M e j and E sn incorporate the range 
inferred from observations of "normal" SNe, but 
do not include hypernovae. 

• We assume that the L x — E relationship used in 
£|3 . 21 is accurate to better than two orders of mag- 
nitude. As a result, in H4.1.U we assume that the 
current spin-down luminosity of the neutron star 
in G328.4+0.2 is one of the following: 

E = 0.017, 0.17, 1.7, 17, and 170 x 10 37 ergs s _1 ,(13) 

and in tH.1.21 assume that 1.7 x 10 35 < E < 1.7 X 
10 39 ergs s _1 . 



These are the initial conditions used in both H4.1.1I and 
fc|4.1.21 The remaining input parameters into the model 
are the initial spin-down luminosity Eq and space ve- 
locity v ns of the neutron star, and the method for deter- 
mining the possible values of these parameters is given in 
i Rl.lland ^4. 1.21 Finally, for all trials discussed in S RI. II 
and i j4.1.2l the model begins at a time t = 0.5 years, with 
At = 0.5 years. 

4.1.1. G328.4+0.2 as a Composite SNR 

In this Section, we evaluate the possibility that 
G328.4+0.2 is a composite SNR. To do this, we first 
assume that the outer edge of the radio emission from 



this source denotes the edge of the SNR, and therefore 
-RG328 = flsnr- As a result, for a given value of E sn , 
M e j, and n, we use the model for the evolution of a SNR 
discussed in Appendix [A] to calculate the current age of 
G328.4+0.2, t now . With this value of t aow and assumed 
values for the current spin-down luminosity of the neu- 
tron star, E, and ro, we are able to calculate both the 
initial spin-down luminosity Eq and initial period Pq of 
the neutron star in G328.4+0.2 using Eqs. (TAT]) and |[X9]) . 
respectively. 

Using this procedure, we ran our model using all pos- 
sible combinations of the input parameters (tq, E, E sn , 
M e j , and n) given in §4.11 for a total of 540 different com- 
binations. To see which combination of these input pa- 
rameters provide a plausible explanation for G328.4+0.2, 
we require the following: 



• Criterion 1: 



< 2000 km s , since a neu- 



tron star with a higher velocity than this is ex- 
tremely implausibl e based on pulsar observations 
dHobbs et al.ll2005h . 

• Criterion 2: Pq > 2 ms, the minimum rota- 
tion period of a young proto-neut ron star before 
it breaks up (jGoussard et al.l lX998). 

• Criterion 3: The PWN is smaller than the SNR 
for all t < t now . While it is possible that PWN 
could expand to fill the entire SNR, it is not con- 
sidered likely, and is contrary to the model as- 
sumption that the PWN does not affect the evolu- 
tion of the SNR. Additionally, we also require that 
flpwn(*now) > 0.67i? snr (i now ), due to the large ob- 
served size of the flat-spectral index radio emission 
which is produced from the PWN, as described in 



• Criterion 4: G328.4+0.2 is currently in the Free- 
Expansion or Sedov- Taylor phase of its expansion, 
i.e. i now < t ra d, where i rac i, the age when a SNR 
goes radiative, is defined in Eq. IA4I Once a SNR 
has entered its Radiative phase, it is expected that 
the radio emission from the SNR be confined to 
thin, brig ht, filaments like those obs erved in SNR 
G6.4-0.1 (|Mavromatakis et all [2004D - which are 
not observed in G328.4+0.2, or that the SNR is 
radio-quiet. 

• Criterion 5: The PWN/RS collision has already 
occurred, as described in i j3.2[ and the PWN has 
been compressed as a result of its collision with the 
RS. This is required to explain the Central Bar, as 



described in 



This requires that n pwn < at 



some point in the past - which can only occur at 
a time t > t co i, the time when the PWN and RS 
collide. 

• Criterion 6: The Central Bar created by the 
compression is still observable. This is satisfied 
if either the PWN is currently being compressed, 
v pvm {tnow) < 0, or if the compression ended suffi- 
ciently recently such that it can be observed. The 
Central Bar is believed to be formed by both a 
pressure and a mag netic field enhancement at the 
center of the PWN (|Bucciantini et al.ll2003l ). As a 
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result its observable lifetime is the synchrotron life- 
time of electrons accelerated by the magnetic field 
enhancement. Therefore, we assume that the life- 
time of the central bar is the synchrotron age of 
the accelerated electrons, T syncn , equal to: 

Ts y „ch = 3 x lO^-^B-V* years, (14) 

where v is the observed frequency, in units of Hz, 
and -Bpwn is the strength of the magnetic field in- 
side of the PWN, in units of G. With this infor- 
mation, in the case that w pW n(inow) > we deter- 
mine if the central bar is still observable by eval- 
uating r sync h at 22 GHz (since this is the highest 
frequency at which the Central Bar is observed; 
iJohnston et al.l l"2004h at the time when the com- 
pression ends, can be derived 
using the minimum energy estimate. This criterion 
is satisfied if T S y ncn ^> t now tic— cxp? where t rc — CX p 
is the time when the compression phase ends. 

That the above criteria fall into two categories: criteria 
required for physical plausibility (Criteria 1-3) and those 
which depend on our interpretation of the radio and X- 
ray properties of G328.4+0.2 (Criteria 4-6). 

Of the 540 possible combinations of the input parame- 
ters, only one passes all six criteria. The predicted SNR, 
PWN, and neutron star properties of this scenario are 
given in Table and the behavior of R pwn as a function 
of time is given in Figure[51 In this scenario, G328.4+0.2 
is quite young, ~ 4900 years old, and the energy in- 
jected into the SNR by the neutron star is similar to 
the kinetic energy of the SN explosion (~ 10 51 ergs). 
However, in this scenario, as shown in Fig. [8l the pre- 
dicted compression is very small; when the compression 
begins, i? pwn = 3.838 pc, and when the PWN begins to 
re-expand into the SNR, R pwn — 3.834 pc. This small 
decrease is not surprising given that, as shown in Table 
El the total energy inputed into the PWN by the pulsar 
(E psr ) is very close to the kinetic energy of the SN ejecta 
(E sn ). This negligible decrease in the volume of the PWN 
is unlikely to form a central bar as prominent as the one 
observed (Fig. [3]), and therefore we feel is unlikely to be 
the correct explanation for G328.4+0.2. 

4.1.2. G3284+0.2 as a PWN 

In order to evaluate if G328.4+0.2 is a PWN, we use 
the model presented in §4] to determine the earliest time 6 
(t now ) at which a PWN powered by a neutron star with 
a given initial period Pq reaches the observed size of 
G328.4+0.2 (i?pwn = 12.5 pc) if it is expanding into as 
yet unseen SNR formed by ejecta with initial mass M c j 
and kinetic energy E sn which exploded in a constant- 
density ambient medium with number density n. To 
consider all reasonable cases, we ran our model using all 
combinations of the values of To, E sn , M e j, and n given 
in £|4. H as well as P = 5, 10, 25, 100 ms for a total of 
432 different trials. In this scenario, since it is not pos- 
sible to determine an independent estimate of the age 
of the system, it is necessary to assume a value of Pq. 
To determine which of these combinations are possible 
explanations for G328.4+0.2, we required that: 

6 Due to oscillations in radius the PWN undergoes after its col- 
lision with the RS, it can reach the current size at multiple times. 



• Criterion 1: u™ in < 2000 km s" 1 . Since, in 
this scenario, we have no prior estimate of the 
age of G328.4+0.2, when we run our model we as- 
sume that v ns — 0. This does not effect our re- 
sults because r ns , as measured in §4.1[ is less than 
i? pwn = i?G328: and therefore the neutron star is 
always injecting energy into the PWN as it does if 
v ns = 0. Once, for a given set of input parameters 
we have determined t now , we calculate u™ m using 
Equation ©. 

• Criterion 2: The current spin-down energy of the 
neutron star in G328.4+0.2 is between 0.017 < 
E,37 < 170, where E >37 = E/10 37 ergs. This is 
based on the work done in §3.2[ and is consistent 
with the initial values of E used in ^4.1.11 Since 
in this scenario we have no estimate of the age of 
G328.4+0.2, we are unable to assume a value for 
E of the central neutron star and then calculate its 
initial spin-down luminosity, as we did in §4.1.11 

• Criterion 3: Rp Wn < i? S nr for all times t < t now , 
as explained in §4.1.11 

• Criterion 4: The PWN has already collided, and 
has been compressed by, the RS, as explained in 

3321 

• Criterion 5: The Central Bar created by the 
compression of the PWN is still observable. The 
method of determining if this is satisfied is the same 
as the one used in ^4.1.11 

• Criterion 6: The PWN must have been able to 
form RT instabilities after the PWN/RS collision. 
As in §4.1.1[ we implement this requirement by re- 
quiring that P pwn > P sm (Rp wn) fo r some t > t co \. 
Additionally, as explained in W3.21 in order for the 
PWN to create Filamentary Structure C it must 
currently be unstable to R-T instabilities - requir- 
ing that P pwn > P sn r(^pwn) nOW. 

• Criterion 7: As explained in §221 the ob- 
served Outer Protrusions in the radio require that 
the PWN currently be expanding into the SNR, 

^pwn(^now) ^* 0. 

• Criterion 8: G328.4+0.2, must have only under- 
gone one compression/re-expansion cycle. As ex- 
plained in ^3.21 numerical simulations of PWN in- 
side SNRs finds th at the PWN is disru pted after 
the first such cycle (Blondi n et al.ll2~00lD . 

It is important to note that, if G328.4+0.2 is a PWN, 
then the radio and X-ray observations provide little in- 
formation on the evolutionary phase of the (unseen) SNR 
and no information on the current ratio of the PWN and 
SNR radii. 

Out of the 432 possible combinations of the input pa- 
rameters, only five satisfy all ten of the above criteria, as 
listed in Table 2) While the neutron star appears to be 
inside the PWN, it is possible that this is just a projec- 
tion effect. To evaluate the possibility that the PWN in 
G328.4+0.2 has already entered the Relic PWN phase of 
its evolution, we calculated u™ m,n , defined as: 
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If v ns > iC n ' n , then the G328.4+0.2 is a Relic PWN, if 
not, then it is still in the Collision with the RS phase of 
its evolution. 

As shown in this Table, the predicted properties of 
G328.4+0.2 vary substantially if G328.4+0.2 is inside a 
Sedov or Radiative SNR. In the Sedov case, G328.4+0.2 
is quite young, and progenitor SN explosion had a normal 
explosion energy but a low ejecta mass (M e j ~ 1-M©), 
and it occurred in a low density environment. Ad- 
ditionally, the neutron star formed in this explosion 
was spinning rapidly, has a low surface magnetic field 
strength (P ns < 10 12 G), and a high space velocity 
(Vna £ 800 km s" 1 ). In the Radiative Case, G328.4+0.2 
is substantially older, and the progenitor SN explosion 
was a low kinetic energy (E sn ~ 3 x 10 50 ergs) and high 
ejecta mass expanding. The neutron star in this case 
was born spinning somewhat slower and has a normal 
magnetic field strength for a young neutron star. 

With the information presented in Table (3 it is possi- 
ble to further refine the expected SNR and PWN prop- 
erties. As argued in £14.1.11 the prominence of the central 
bar argues that, during the compression stage, the vol- 
ume of the PWN decreased significantly. Though it is not 
possible at this time to quantify the compression needed, 
an examination of Table [5] shows that for only two mod- 
els, ST 2 and Rad 1, did the volume of the PWN decrease 
by more than 10% - and therefore these two models are 
the most probable descriptions of G328.4+0.2. In the 
case of Rad 1, «™ n ' n ~ 100 km s _1 is significantly less 
then the average neutron star velocity, v ~ 400 km s _1 
(jFaucher-Giguere fc Kas"pll2006t iHobbs et al.ll2005l ). im- 
plying that the PWN is in the Relic PWN phase of its 
evolution. Since the sound speed inside a Radiative SNR 
is quite low, ~ 100 km s _1 , we expect that the PWN in 
the Rad 1 scenario would have a bow-shock morphology. 
Since there is no clear evidence for this in the X-ray or 
radio emission from G328.4+0.2, this suggests that ST 2 
is a better fit to the data. 

The conclusion that ST 2 is an accurate description 
of G328.4+0.2 is supported by circumstantial evidence 
as well. For this scenario, the expected radius of the 
ter mination shock ar ound the neutron star, r ts , defined 
as (jSlane et al.ll2004h : 

ri = sJ— - (16) 

^' lu pwn 

assuming a spherical wind, is r ts ~ 0.6 pc, which corre- 
sponds to an angle of 6> ts ~ 8^17" - a distance which is 
comparable to the offset between Clump 1 and Clump 2 
derived in £12.1.11 Another interesting feature for this 
model is that, as shown in Fig. [4) the radius of the PWN 
at the time of re-expansion is similar to that of the outer 
parts of the central filamentary structures discussed in 
£12.21 While this correlation might be coincidental, this 
could imply that the hydrodynamic instabilities formed 
at the PWN/SNR interface during the re-expansion dis- 
rupted the shell of material swe pt up by the PWN - 
consistent with the simulation of iBlondin et al.l {2001). 
While not definitive, these two pieces of evidence argue 
that ST 2 is a reasonable description of G328.4+0.2. 

The properties of ST 2 are given in Table El and the 
evolution of i? pwn is shown in Fig. [9] It is interesting to 
note that in this scenario, the PWN collides with the RS 
at a time t co \ < tq (t co \ s=s 850 years) so energy injection 



by the neutron star into the PWN after the PWN/RS 
collision is important to the PWN's evolution during this 
stage. It is important to note that this model predicts 
that the neutron star powering G328.4+0.2 is the most 
energetic neutron star in the Milky Way, as well as one 
of the fastest. Given that G328.4+0.2 is largest and has 
the highest radio luminosity of any known PWN, it is not 
surprising that it was formed by such a powerful neutron 
star. Finally, the age and E pr edicted by this metho d 
are similar to those predicted bv lGaensler et al.l (2000). 

In order to better understand the limitations of the 
approach in determining the properties of the neutron 
star and SNR in G328.4+0.2, we have run the model 
presented in £j4]over a finer grid of parameters and eval- 
uated the resulting PWN evolution using the same cri- 
teria as above. In Fig. [TD1 we show which values of Po 
and B ns pass all of this criteria for three different kind 
of SN explosions: P sn = 10 51 ergs and M cj = 1 M Q , 
E sn = 3 x 10 51 ergs and M ej = 1 M , and E sn = 

4 x 10 51 ergs and M e j = 3.25 M Q , assuming an ambi- 
ent density with n = 0.03. The first set of SN parame- 
ters corresponds to ST 1, the second to ST 2, and third 
to a higher ejecta mass SN explosion is compatible with 
a neutron star with the same parameters as ST 2. For 
the first case, we find that a wide range of Po values 
are allowed but that B ns < 10 12 G. In fact, for this set 
of SN parameters a Po ~ 10 ms, B ns ~ 8 x 10 11 G neu- 
tron star results in a PWN which is compressed a similar 
amount as in ST 2. In the second case, we find that Po 
and P ns are tightly constrained around Pq sa 5 ms and 
B DS » 5 x 10 11 G. In the third set of SN parameters, we 
find that Po < 6 ms, but that P ns spans a wide range of 
values, - 10 11 - 2 x 10 12 G. 

To determine the allowed values of P sn and M e j, we 
followed the same procedure as above using two different 
sets of neutron star parameters: Po = 5 ms and P ns = 

5 x 10 11 G (the neutron star parameters in the ST 1 and 
ST 2 scenarios), and P = 10 ms and P ns = 8 x 10 11 G. 
For the first case, only models with P sn ~ 1— 4x 10 51 ergs 
and M e j ~ 0.5— 3.5M Q satisfy the criteria above - though 
a substantial compression of the PWN requires P sn > 
2 x 10 51 ergs. In the second case, we find that only 
models with E sn < 10 51 ergs and M C j ~ 0.5 — 3.5M Q are 
allowed. 

While this error analysis shows that the method used 
above to determine the properties of the neutron star 
and SN explosion which formed G328.4+0.2 is unable to 
do so to much better than an order of magnitude, the 
different combinations values of Po, P n s, Psn, and M j 
which are allowed predict different physical properties for 
G328.4+0.2 which are testable with further observations. 
For example, in the case of Po = 10 ms, P ns = 8 x 10 11 G, 
E sn = 1 x 10 51 ergs, and M ej = 1 M , G328.4+0.2 is 
~ 13, 000 years old, twice the age predicted in the ST 2 
model as shown in Table O and as a result the required 
velocity of the neutron star is significantly lower, i>™ m ~ 
500 km s _1 . The predicted period for the neutron star 
in this scenario is also significantly slower than required 
by the ST 2 scenario, P ~ 24 ms, with a value of E 
approximately an order of magnitude lower than that 
in the ST 2 scenario. The termination shock radius for 
this set of parameters is ~ 5", considerable smaller than 
the ~ 8" for the ST 2 scenario and detectable with the 
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Chandra X-ray Observatory. 

5. CONCLUSIONS 

In this paper, we first presented new X-ray ( §2.1|) and 
radio (§ §2.21 I2.3|) and observations of Galactic non- 
thermal radio and X-ray source, G328.4+0.2, from which 
we infer the current properties and evolutionary history 
of this source f ^3.2p . We then presented a simple hydro- 
dynamic model for the evolution of a PWN inside a SNR 
(§4j, which is used to determine which values of E sn , 
M e j, n, Pq, and B ns are able to reproduce the properties 
discuss ed in §321 if G328.4 +0.2 was a Composite SNR 
( gXP or a PWN ( H4.1.2|) . As a result of this analysis, 
we determine the G328.4+0.2 is a PWN inside an unde- 
tected SNR. Though we are not able to precisely deter- 
mine the properties of the SN explosion and the neutron 
star which have created this system, our analysis implies 
that the neutron star in G328.4+0.2 was born with an ini- 
tial period Pq < 10 ms, has a lower than average surface 
dipole magnetic field strength, and has a higher than av- 
erage spatial velocity v ns > 400 km s _1 . We assume de- 
termining that the SN explosion which created the neu- 
tron star had a normal explosion energy, E sn ~ 10 51 ergs, 
but a relative low ejecta mass, M e j < AMq. Future X- 
ray and radio observations can significantly decrease this 
uncertainty, particularly if they are able to either detect 
pulsations from the neutron star or continuum X-ray or 
radio emission from the currently undetected SNR in this 
system. 

While we are not able to definitely determine the ini- 
tial period (Pq) or surface magnetic field strength (B ns ) 
of the neutron star, nor the kinetic energy (E sn ) or ejecta 
mass (M e j ) of the progenitor SN explosion, the estimates 
quoted above are of interest. Our non-detection of the 
pulsar via radio pulsations is not particularly constrain- 
ing, due to the very large distance of the PWN. The low 



magnetic field but rapid initial period predicted for the 
neutron star in G328.4+0.2 has implications for mod- 
els concerning the origin of neutron star magnetic fields. 
For example, according to t he a — Q dynamo model of 
([Thompson fc D uncan 1993), neutron star born spinning 
rapidly (Pq < 5 ms) should have a strong dipole compo- 
nent to their surface magnetic fields (B ns ^> 10 12 G). 
If the ST 2 scenario proves to be correct, then the low 
magnetic field of the neutron star in this system would 
be a problem for such a model. Additionally, the low 
ejecta mass inferred in this scenario requires that the 
progenitor of this system was either a single, massive 
star (M > 35Afc) w hich exploded in a Type Ib/c SN 
(|Wooslev et al.|[l995f) . or was initially in a binary sys- 
tem. 

Finally, the method used in this paper to study 
G32 8.4+0.2 is comp lementary to other methods used 
(e.g. IChevalierll2005f ) to infer the initial period and mag- 
netic field strength of other neutron stars in young PWN 
as well as the properties of the SN explosion in which 
they were formed, and is easily applicable to other such 
systems. 
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Table 1 

Spatial components of the X-ray emission from G328.4+0.2 



Component 


Parameter 


Value 


background 


constant 


09 +u - Ud ' 

u - uy -0.07 


Clump 1 


ro 


,». 

l'/i+V 

-07 7 




Position 15 h 55 m 26?6 
.4 
a 


3+ o=05' -53°18'02"7+ ';', 6 
°-9-o°3 



10" 

+0?2 _5 3 o 17 / 54 // 6 +3'.'4 



Clump 2 



FWHM 
Position 

e 



.4 



19" 



15 h 55 m 27!5 



-ori' 
0.4 

280 

0.8 



+0.2 
0.3 
o+20° 

80° 
+0.3 
0.2 



Diffuse 



FWHM 
Position 

e 



67 



-11 



15 h 55 m 26!6 +0 s 2 , -53°17'48'.'3 +5 „ 

„+0.1 

-0.1 
+ 10° 



0.4 
130° 
0.4 



10" 
+0.1 
0.1 



Note. - Results from the spatial fit to the X-ray emission from G328.4+0.2, as described in £|2.1,ll Fitting was done using the Sherpa 

f 2 \~ a 

software package, and the errors reflect the 90% confidence level. Clump 1 was fitted to a 2D Lorentzian, defined as f(r) = A I 1 + ^ J , 

where f(r) is the expected number of counts at radius r away from the center of the source, A is given in counts, and the core-radius ro is in 
arc-seconds. The background component is given in counts pixel -1 . Both Clump 2 and Diffuse were modeled with elliptical 2D Gaussians, 
where the full width, half maximum (FWHM) is given in arc-seconds, the ellipticity is e, defined as 1 — b/a, where b and a are, respectively, 
the major and minor axis of the source, the position angle 9 is given in degrees counterclockwise from north, and the amplitude A is given 
in counts. 



Table 2 

Spectral Fits to Mosl + Mos2 + pn Data 



Parameter 


Value 


Model 


phabs * pow phabs * bbodyrad phabs * bremss phabs * ray 


T or kT 
Absorbed Flux 
Unabsorbed Flux 


11.6+2 x 1q22 6 - 4+ ! x 1q22 10.4+^ x 10 22 7.8+i x 10 22 
2.0«-j 1.6±g;i 9.5±§ 62.9i 3 >r 

4.6+3 x 10-13 4 - 4+ o 5 x 10 ~ 13 4 - 5+ °'5 x 10-13 4 - 9+ 2'e x 10-13 
1.9±g;| x lO" 12 6.5±?;! x lO" 13 1.2±g;J x 10" 12 9.6til x lO" 13 


x 2 

Reduced x 2 


103.3/131 98.4/131 101.7/131 122.8/131 
0.81/128 0.77/128 0.79/128 0.96/128 



Note — Results from joint fits to the MOSl, MOS2, and PN spectrum between 0.5 and 10 keV. The model phabs * pow refers to a power-law 
attenuated by interstellar absorption, phabs * bbodyrad refers to a blackbody attenuated by interstellar absorption, phabs * bremss 
refers to a Bremsstrahlung source spectrum attenuated by interstellar absorption, and phabs * ray refers to a Raymond-Smith thermal 
plasma spectrum attenuated by interstellar absorption. In the table, Nj^ is given in cm -2 , kT in keV, and flux in ergs cm -2 s — 1 , both the 
absorbed and unabsorbed flux are calculated between 0.5 and 10 keV, and the errors represent the 90% confidence level. 
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Table 3 

Expected Properties of G328.4+0.2 if it is a Composite SNR 



Supernova Remnant Properties Pulsar Wind Nebula Properties Neutron Star Properties 



Parameter 


Value 


Parameter 


Value 


Parameter. 


Value 


E sn 


1 X 10 51 ergs 




4.8 X 10 50 ergs 


E 


2.5 X 10 40 ergs s" 1 


Mej 


1 M 


^pwn 


8.7 parsecs 


Po 


3.8 ms 


Phase 


Sodov- Taylor 


M pwn 


5.3 M e 


TO 


1730 years 


Psnr (^ipwn ) 


1.4 X 10~ 9 dynes 


^pwn 


2.0 X 10~ 9 erg cm -4 


B ns 


1.1 x 10 12 G 


^cj (^pwn ) 


400 km s" 1 


"pwn 


700 km s" 1 


7 ,min 
ns 


1300 km s" 1 


£now 


4900 years 


rta 


0.5 parsecs 


E 


1.7 X 10 39 ergs s _1 


n 


0.32 cm 3 


Ots 


6" 7 


P 


7.5 ms 










P 


1.8 x 10~ 14 s/s 










Tc 


6700 years 












1.0 X 10 51 ergs 



Note. - The values in bold are model assumptions, while the others are predicted by the model presented in [|4] P S nr(fipwn) is the 
pressure inside of the SNR just outside of the PWN, v c j (R pvln ) is the velocity of material inside the SNR just outside of the PWN, E pW n is 
the internal energy of the PWN, R pvin is the radius of the PWN, Mj^" 1 is the mass of material swept up by the PWN, P pW n is the internal 
pressure of the PWN, u pW n is the exp ansi on velocity of the PWN, rt B is the radius of the termination shock around the neutron star inside 
the PWN, calculated using Equation 1161 6ts is the angular size of this feature assuming d\i = 1, B nB is the dipole magnetic field of the 
neutron star in G328.4+0.2 according to this model, P is the period of the neutron star in G328.4+0.2, P is the period-derivative of the 
neutron star, t c is the characteristic age of the neutron star, defined as r c = P/(2P), and -E pS r is the total amount of energy injected by 
the neutron star into G328.4+0.2 for t < t novJ . All values are given for t = i now unless otherwise noted. 



Table 4 

Scenarios for G328 as a PWN inside an Undetected SNR 



Scenario # 


TO 


Po 


-Esn, 51 


M cj 


n 


tnow 


SNR Phase 


Rsnr 


B n s,12 


min 
u ns 


min.il 


ST 1 


1730 


5.0 


1.00 


1 


0.03 


5100 


Sedov- Taylor 


19.8 


0.5 


1300 


2400 


ST 2 


1730 


5.0 


3.16 


1 


0.03 


6500 


Sedov- Taylor 


28.0 


0.5 


1000 


1900 


Rad 1 


430 


10.0 


0.32 


10 


0.32 


84200 


Radiative 


28.4 


2.0 


100 


100 


Rad 2 


770 


10.0 


0.32 


10 


0.32 


52400 


Radiative 


24.8 


1.5 


100 


200 


Rad 3 


770 


10.0 


0.32 


10 


1.00 


101400 


Radiative 


22.1 


1.5 


100 


100 



Note. - Values for to, Po, -Ban, M e j, and n that satisfy the criteria listed in £14.1,21 for G328 being a PWN inside an undetected SNR. 
The value of to is given in years, Po in ms, Esn, 51 = E sn /10 51 ergs, M e j in solar masses, n in cm , inow in years, Rsnr in parsecs, 
Bns = Bns, 12 X 10 12 G, v^ n is in km s -1 , and l%g in ' is also in km s _1 . 



Table 5 

Compression/Expansion properties of G328.4+0.2 if it is a PWN 



Scenario # 


t(u P wn = 0) 


Central Bar Lifetime 


Rpwn (f pwn — 0) 




ST 1 
ST 2 


1648.5, 2432.0 
969.0, 1837.0 


63930 
35250 


8.20, 7.96 
8.36, 6.34 


0.91 
0.44 


Rad 1 
Rad 2 
Rad 3 


12772.0, 26541.0 
13336.5, 21112.5 
9761.5, 11035.0 


243480 
256230 
102640 


8.53, 7.55 
8.48, 8.27 
5.72, 5.72 


0.69 
0.93 
1.00 



Note. - In this table, values Scenario # correspond to the values of to, Po, E Bn , M e j, and n given in Table|4] 4(i> pW n = 0) and the Central 

V-re-cxp 

Bar lifetime are given in years, and R pW n(f P wn = 0) is given in parsecs. {SiSpivai is the ratio of the volume of the PWN at re-expansion 

'pwn 

and compression. 
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Table 6 

Expected Properties of G328.4+0.2 if it is a PWN (ST 2 scenario in Table HJ) 



Supernova Remnant Properties Pulsar Wind Nebula Properties Neutron Star Properties 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 


E sn 


3.2 X 10 51 ergs 


^pwn 


3.2 X 10 50 ergs 


Eo 


1.5 X 10 40 ergs s" 1 




1 M 


^pwn 


12.5 parsecs 


Po 


5 ms 


Phase 


Scdov- Taylor 


, ,pwn 
Jw sw 


0.4 M Q 


TO 


1730 years 


PsnrC^pwn ) 


3.6 x lCT 10 dynes 


P 

1 pwn 


4.4 x 10~ 10 dynes 


B n s 


5 x 10 11 G 


^cj (^pwa ) 


460 km s" 1 


Wpwn 


790 km s" 1 


7 .min 

ns 


990 km s" 1 


Age (Ww) 


6500 years 


rts 


0.6 parsecs 


E 


6.4 X 10 38 ergs s _1 


n 


0.03 cm 3 


Sts 


7'.' 7 


P 


10.9 ms 










P 


2.1 X 10~ 14 s/s 










Tc 


8200 years 










Epsr 


6.3 X 10 50 ergs 



Note. - The model assumptions are given in bold. P S nr (-Rpwn ) is the pressure inside just outside of the PWN, f c j(iJ P wn) is the velocity 
of material inside just outside of the PWN, E pvln is the internal energy of the PWN, i? pW n is the radius of the PWN, A/J^" 1 is the mass of 
material swept up by the PWN, P pW n is the internal pressure of the PWN, v pvn is the expansion velocity of the PWN, rts is the radius of 
the termination shock around the neutron star inside the PWN, 6t B is the predicated angular radius of this feature assuming di7 = 1, S ns 
is the predicted dipole magnetic field of the neutron star in G328.4+0.2, P is the predicted period of the neutron star in G328.4+0.2, P 
is the predicted period-derivative of the neutron star, r c is characteristic age of the neutron star, and E pBI is the total amount of energy 
injected by the neutron star into G328.4+0.2 for t < t n ow- All values are given for t = t n ow unless otherwise noted. 
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-17' 



Background Region for Spectral Analysis 



Background 
Point Source 



Clump 2 Diffuse 




Clump 1 



Source Region for Spectral Analysis 



o 
o 
o 



_L 



15 h 55 m 35 s 30 s 25 s 20 s 

Right Ascension (J2000) 



15 s 




38 s 26 s 
Right Ascension (J2000) 



Fig. 1. — Top: Exposure normalized, vignette corrected MOSl + MOS2 image of G328.4+0.2 smoothed by a 5" Gaussian. The white 
contours indicate 20, 40, 50, 70, and 90% of the peak X-r ay flu x in the smoothed image, while the boxes indicate the background and 
source regions used for the spectral analysis described in £12.1,21 The background point source labeled in t he im age was excluded from 
the background region. Additionally, the labels in this plot point to the morphological features discussed in £|2,l.ll Bottom: Unsmoothcd 
nor malize d, vignette corrected MOSl and mos2 image of G328.4+0.2 overlaid with the regions used for the Hardness Ratio analysis discussed 
in fl2~l~2l 
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Fig. 2. — Mosl, MOS2, and PN spectrum of G328.4+0.2 overlaid with the absorbed power-law model whose parameters are given in 
Table [2] 
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Fig. 3. — 1.4 GHz image of G328.4+0.2, overlaid with X-ray contours in green which represent 20%, 35%, 90% of the peak flux in the 
smooth X-ray image shown in Fig. [T] The beam size of this image is 7'.'0x5'.'8, and is sho wn in the lower left-hand corner of the image. 
The labels indicate examples of the different radio morphological features discussed in £12.21 



20 



-53°14' 



- 16' 



-ia 1 - 



"i j 1 1 s r~ 



"! f 1 f~ 



~t 1 1 t- 




J_ 



15 h 55 m 50' 



_L 



_L 



15 



e 

10 s 
m 

S 



40* 30 s 

Right Ascension (J2000) 



20 s 



10 s 



Fig. 4. — 20cm radio i mage of G328.4+0.2 (same data as shown in Fig. [3), with a color scale chosen to enhance the visibility of Filamentary 
Structure B discussed in i|2.2l The yellow circle indicates the size of PWN predicted in the ST 2 model listed in Table [5] when it re-expanded 
after the initial compres sion by the SNR rever se shock, and the white cross indicates the center of G328.4+0.2 (15 h 55 m 33 s , -53°17'00"; 
J2000) as determined bv lGaensler et al.l l l2000h . 
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15 h 55 ra 50%0 , 30*20* 40* 30 s 20 s 40 s 30 s 20 s 40 s 30 s 20 s 



HA. KA KA. RA 

Fig. 5. — Spectral tomography images of G328.4+0.2, as described in i|2.2,ll The spectral index a is given in the upper left hand corner 
of each image, where S v oc u a . 



SNR Forward Shock 




Fig. 6. — Diagram of a Composite SNR in the Free Expansion stage of its evolution. In this image, the ratio between the thickness 
of the mass shell surrounding the PWN and the radius of the PWN is 1/24, as determined by Ivan der Swaluw et al.l (t200lT) ■ a nd the 
radio of the SNR Forward Shock, Contact Discontinuity, and Reverse Shock radii are equal to the values given in lCheva licr ( 1982) for his 
n = 9, s = case. The colors denote the nature of the material within each region. 
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Fig. 7. — -Rpwn/flsnr for Models A (solid), B (long-dashed line), & C (short-dashed line) in Blondin ct al. (2001). The top plot shows 
the result of the model presented in [|4] while the bottom is a reproduction of Fig. 3 by Blondin e t al. (2001), reproduced by permission of 
the A AS. 
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35 
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SNR Radius 
Reverse Shock Radius 
Neutron Star position 
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0.1 




100 



1000 
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Fig. 8. — The radius of the PWN, SNR, and SNR reverse shock as well as the location of the neutron star as a function of time if 
G328.4+0.2 is a composite SNR. The vertical line indicates the current age of the system, and the properties of this system are given in 
Table [3] 



PWN Radius 
SNR Radius 
Reverse Shock Radius 
Neutron Star position 




Fig. 9. — The radius of the PWN, SNR, and SNR reverse shock as well as the location of the neutron star as a function of time for the 
favored (ST 2; Table [4} scenario if G328.4+0.2 is a PWN. The vertical line indicates the current age of the system. 
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Fig. 10. — The results for varying Pq and B ns for a E sn = 10 51 ergs, M oj = 1M Q (top), E sn = 3 X 10 51 ergs, M cj = IMq (middle), and 
E S n = 4 X 10 51 ergs, M e j = 3.25 M(T) (bottom) SN explosion, assuming n = 0.03 cm -3 . The small black squares indicate models which failed 
the criteria described in £|4.1,2I while the colored circles indicate scenarios which passed. The color represents the Compression Fraction 
of the PWN, defined as the ratio of the PWN's volume at the beginning and end of the compression stage. The Compression Fraction of 
the ST 2 case given in Table [4] is 0.44 (light blue on this color scale), and lower values correspond to a more substantial compression. The 
star indicates the position of a Po = 5 ms, B nB = 5 X 10 11 G neutr on s tar, while the large square indicates the position of a Pg = 10 ms, 
Bns = 8 X 10 11 G neutron star - the two neutron stars used in Fig. 1111 
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Fig. 11.— The results of varying Esn stud -^^cj for 9> -Po — ^ ms, — 5 x 10 11 G (top) and P = 10 ms, B nB = 8 x 10 11 G (bottom) 
neutron star, assuming n = 0.03 cm -3 . The black squares indicate which scenarios failed the criteria described in £14.1,21 while the colored 
circles indicate those that passed, with the color representing the Compression Fraction (defined as the ratio of the PWN's volume at the 
beginning and end of the compression stage) of the PWN. The star indicates a E sn = 1 X 10 51 ergs, M e j = 1 Mq SN explosion, the square 
indicates a E sn = 3 X 10 51 ergs, M e j = 1 Mm SN explosion, and the circle indicates a E an = 4 X 10 51 ergs, M e j = 3.25 Mq SN explosion - 
the SN explosion parameters used in Fig. 1101 
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APPENDIX 

EQUATIONS FOR THE HD MODEL FOR THE EVOLUTION OF A PWN INSIDE A SNR 

In this Appendix, we provide many of the details concerning the properties of the neutron star, PWN, and SNR 
needed to simulate the HD model for the evolution of a PWN inside a SNR described in 

For the SNR, we assume that the initial ejecta density profile consists of a constant density core surrounded by a 
p oc r ~ 9 envelope - the standard assumption for a SNR produced by a Type-II SNR (jBlondin et alj|200ll : IChevalierl 
1982) - and that the ejecta is expanding ballistically (u e j = r ? -Jt). The boundary between the constant density core 
and the outer ejecta envelope has a velocity v core , defined as (jBlondin et al.|[200ll i: 



20 E^ 1/2 
9 M cj 



(Al) 



where E sn is the explosion energy of the SN and M e j is the ejecta mass. As a result, the density p COIC of the ejecta 
core is (jBlondin et al.ll2001h : 

Pco re (*) = ^ E sn u- 5 re r\ (A2) 

As the SNR expands, it sweeps up and shocks the surrounding interstellar medium (ISM). This swept-up material has 
a higher pressure than the cold ejecta driving the expansion of the SNR, and as a result drives a reverse shock (RS) 
into the SN ejecta. In between the outer edge of the SNR, which marks the location of the forward shock (FS) and 
the RS is a contact discontinuity which separates the shocked ISM from the ejecta shocked by the RS. The pressure 
inside the SNR at r < r rs , where r rs is the radius of the RS, is assumed to be zero. A diagram of this is shown in 
Fig. [3 Since we assume that both the SN ejecta and the shocked ISM behave as a 7 = 5/3 perfect gas, the sound 
speed c s of this material is: 




(A3) 

When the RS is still in the ejecta envelope, we determine the pressure, velocity , and density profiles on the material 
between the RS and FS using the self-similar equations given by IChevalierl (|1982f ). evaluating them for the n = 9, s = 
case. How ever, when the R S enters the constant de nsity ejecta core, it i s no l onger possibly to apply the self-similar 
solution of lChevalierl (119821), and we use the work of iTruelove k, McKed (j!999l) to determine the radius of the RS and 
the results given bv lBandieral |l984) to determine the pressure, velocity, and density profil es between the RS a nd FS . 
It is also necessary to model the radius of the FS (i? snr ), which we do using the work of ITruelove fc McKed (1999). 
This is valid while the SNR is in the Free Expansion and Sedov- Taylor p hases of its evolutions. After the SNR goes 
radiative, which occurs at a time t = i ra d defined as (B londin et al.lll998f ): 

i rad «2.9£ 5 f n~* x 10 4 yr. (A4) 

After this point, i? snr oc t 2 ^ 7 . An analytic model for the pressure, velocity, and density distribution of a SNR in this 
phase does not currently exist, and therefore it is difficult to extend our model t o this phase, though if one assumes 
that the interior of the SNR evolves adiabatically, R pwn /Rsnr oc i 075 for t > i ra d (jBlondin et al.|[2~001h . 

For the PWN, as mentioned in SJU we assume that it is a bubble filled with a 7 = 4/3 perfect gas. As a result, the 
internal pressure of the PWN, P pwn is equal to: 

E n 



Ppwn — -zrp , (A5) 

"J •'pwn 

where -Epwn is the internal energy of the PWN and V pwn is the volume of the PWN, defined as: 

Vpwn = g^^pwn (A6) 

where i? pwn is the radius of the PWN. The internal energy of the PWN is determined by the rate of energy injected 
into the PWN by the neutron star (E), and energy loss due to its expansion inside the SNR (E^ n ). For E, we use 
the standard assumption that it is equal to: 

E = E U + -j (A7) 

where t is the age of the neutron star, p is the pulsar braking index (p = 3 for a magnetic dipole), r is the characteristic 
timescale of pulsar spin-down, and Eq is the initial spin-down energy of the neutron star. Both m and Ep dep end on 
the physical properties of the neutron star, with tq denned as (jBlondin et al.|[200lt IShapiro fc Teukolsky|[l983[) : 



3c 3 IP 2 

T ° = Ai, 2 B 2 mlm 2 a {M) 
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where I is the neutron star's moment of inertia, Pq is the initial spin period, B ns is the magnetic field of the neutron 
star, R ns is the radius of the neutron star, a is the angle between the neutron star's rotation axis and magnetic field, 
and (jBlondin et al.ll2001h : 



2tt\ 2 1 



E ° = I 7H —f TV ( A9 ) 

\PqJ Mp-i) 



Since the PWN expands adiabatically, -E pwn is equal to 



£$n = -^P- (MO) 



As a result, the change in the internal energy of the PWN over time (E pwn ) is equal to 



E pwn = f^+E [l + -j (All) 

assuming that p — 3. This equation can be solved analytical, and we result that E pwn (t) can be expressed as: 

^.^(MLiiM-^). (A12) 

When we run our model, we use Eq. (|A12|) to determine the initial value of E pvin , but determine E pwn at later times 
using the procedure described in Step 2 in S|4j 

During its free-expansion, the PWN is moving faster than its surroundings, and the mass of the shell surrounding 
the PWN (Af SWiPwn ) is simply: 

fJipw» 

M sw , pwn (i)= / iTrR 2 p ci (r,t)dr (A13) 



where p e j(r) is the density profile of the SNR. After the collision with the reverse shock, if the PWN is moving 
faster than its surroundings we determine the mass of the ejecta shell recently swept up by the PWN and add it 
to the value of M SW)Pwn calculated at the time of the reverse shock collision. If the PWN is moving slower than its 
surroundings, we assume that M SWjPwn remains constant, even if the PWN is being compressed by the surrounding 
SNR. Due to the difference in pressure between the PWN interior to the mass shell and the SNR exterior to the mass 
shell (P sn r( r = -Rpwn))) the mass shell is subject to a force Fapi defined as: 

F AP =47ri?2 wn [p pwn - P snl (R pvjn )}. (A14) 

In this notation, F&p > means that the PWN interior has a higher pressure than inside the surrounding SNR. If the 
PWN has not yet encountered the RS, we assume that P snr (r — R pwn ) — 0. If the mass shell is moving faster than the 
sound speed of the surround ing material (v pvm > c s (i? pwn )), which is the case before the PWN interacts with the RS 
(|Chevalier fc Fransson| [l992l, the mass shell is decelerated by ram pressure, and the total force on the mass swept-up 
by the PWN, F pwn , is: 

F pwn =F A p - 47ri? pwn/ 9 e j(i? pwn )[?; pwn - v cj (R pwn )} 2 . (A15) 

If fpwn < c s: then F pwn = Fap. For t <C to, analytical solutions to these equations give i? pW n °c i 6 / 5 if the PWN is 
still inside the central constant-density core - a result which is reproduced by our numerical implementation of the 
model described in gl 
In this framework, the period P of a neutron star evolves as: 

t \ p+1 

F = Fv\l + 

the period-derivative P evolves as: 

P = 

and the surface magnetic field B ns of the neutron star, assuming p = 3, is 

^. = 1.6 | ^5 mS 1 xH) 9 G (A18) 

where So,37 = Eq/10 37 ergs s _1 , Po.ms is the initial period in ms, and R14 is the radius of the neutron star i?/14 km. 
Additionally, for p = 3 tq is equal to: 

T o = 17.3^|F years (A19) 

^12^14 

where I = lO 45 !^ g cm 2 and the angle between the spin and magnetic field axes of the neutron star is a = 45°. 



p = - p o(l + -) ■ (A16) 



" p+i 



p=£t 1 + ^) (Ai7) 



